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Abstract 

When the vacuum is partitioned by material boundaries with arbi- 
trary shape, one can define the zero-point energy and the free energy 
of the electromagnetic waves in it: this can be done, independently of 
the nature of the boundaries, in the limit that they become perfect 
conductors, provided their curvature is finite. The first examples we 
consider are Casimir's original configuration of parallel plates, and the 
experimental situation of a sphere in front of a plate. For arbitrary ge- 
ometries, we give an explicit expression for the zero-point energy and 
the free energy in terms of an integral kernel acting on the boundaries; 
it can be expanded in a convergent series interpreted as a succession 
of an even number of scatterings of a wave. The quantum and thermal 
fluctuations of vacuum then appear as a purely geometric property. 
The Casimir effect thus defined exists only owing to the electromag- 
netic nature of the field. It does not exist for thin foils with sharp 
folds, but Casimir forces between solid wedges are finite. We work 
out various applications: low temperature, high temperature where 
wrinkling constraints appear, stability of a plane foil, transfer of en- 
ergy from one side of a curved boundary to the other, forces between 
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distant conductors, special shapes (parallel plates, sphere, cylinder, 
honeycomb). 

1 Introduction 
1.1 A brief history 

According to Planck (1900) |T, the energy of a stationary mode of the elec- 
tromagnetic field (with frequency ti;/27r) is quantized with a contribution hw 
per photon. Because of the quadratic structure of the electromagnetic energy, 
a canonical transformation describes each mode as an harmonic oscillator. 
At the classical level, this representation was already known to Rayleigh, and 
Planck guessed rightly at the quantization of the oscillator energy levels. It 
was later discovered that their complete form reads 

eN = nw{N+^), NeN, (1) 

with the number of photons, and with a non-vanishing zero-point energy, 
^hw, which reflects Heisenberg's uncertainty principle. One might have be- 
lieved that the reference vacuum energy would be unobservable. Hendrik 
B G Casimir showed to the contrary in his famous 1948 article i2j that a 
physical force can be generated by vacuum fluctuations (see also jS] , and for 
an historical account, |1]). 

In its original form, the Casimir effect describes an attraction between 
two plane, parallel, perfectly conducting plates, which is explained by this 
occurence of a virtual electromagnetic field in the vacuum separating the 
plates and in the vacuum lying outside them. Indeed, even when no real 
photon is present, the zero-point motion of each mode m of the field yields 
a contribution ^huJm{L) to the energy, which depends on the distance L 
between the plates. The variation with L of the overall zero-point energy 
manifests itself as the Casimir force. It is remarkable that, in spite of the di- 
vergence which appears when summing the zero-point energy over all modes, 
one finds a finite value for this force. Let us remark that a dimensional argu- 
ment allows one to anticipate the right form of Casimir 's result. A quantum 
pressure between two large plates should indeed depend on h, c and on the 
distance L between them. To get pressure units P oc N m~^ = J m~^, with 
/i oc J s, c oc m s~\ and L oc m, one has first to multiply h by c, so as to 
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eliminate time, and then divide by L ^. Hence Pcasimir oc — , and the real 
surprise is the finiteness of the numerical coefficient! 

In recent experiments conducted in 1997-98 ([3 El); the Casimir force 
between a metallized sphere and a metallized plate was finally measured 
definitively with the help of atomic force microscopes. The field is nowadays 
the subject of intense research activity [71 |H]. Measurements of attractive 
Casimir forces are performed for the sphere- plane geometry [HI, for microme- 
chanical torsional devices ^Ul, and for the original Casimir's planar geom- 
etry pi]. Lateral Casimir forces are measured between corrugated surfaces 
[Tallin] (see also [H]). Roughness effects are experimentally important f [T3] 
and references therein). The experimental detection of temperature effects 
is an open problem ^H] • 

1.2 Statement of the problem 

It is somewhat puzzling to regard the Casimir effect as a property of vac- 
uum containing a virtual electromagnetic field. It may look more natural 
to attribute it to the matter of the plates. Actually it is possible through 
Maxwell's equations to express the electromagnetic field, and hence its en- 
ergy, in terms of the charge and current densities of the particles which 
can move within the conducting plates and which are the sources for the 
field. The zero-point motion of these particles then yields a non-vanishing 
interaction energy although the expectation values of the charge and current 
densities vanish. The force between the conducting plates at zero tempera- 
ture can then be interpreted as a result of the interaction between the virtual 
zero-point currents that must exist in the ground state of matter owing to 
Heisenberg's inequality. However, this viewpoint is plagued by the fact that 
the interaction between charged particles is not instantaneous, but retarded. 
On the other hand, the evaluation of the force as a result of this interaction 
would rely on the specific structure of matter. In contrast, Casimir's view- 
point, which focuses on the field rather than on the matter of the plates, 
shows that the effect is nearly independent of the properties of matter (pro- 
vided the plates are good conductors); moreover, as function of the electro- 
magnetic field E, B the energy of this field at a given time is simply expressed 
as 
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in terms of the field at the same time, whereas E(r) and B(r) depend on the 
charges and currents at earher times. 

We wish to study the energy Ef of the electromagnetic field in empty 
regions of space limited by boundaries with arbitrary shape, under circum- 
stances when the expectation value (E), (B) of the quantum field vanishes 
at any point. The energy (j21) may be non-zero for two reasons. 

On the one hand, the uncertainty relations prevent the quantum fluctu- 
ations of E and B from vanishing, since these operators do not commute; 
hence the expression Q for the energy has a positive minimum, the zero- 
point energy of the field. (We shall see that this value is not only strictly 
positive, but in fact infinite; however its variations are finite and physically 
meaningful.) 

On the other hand, the vacuum where the field is considered is bounded 
by walls with which the field can be in thermodynamic equilibrium, at some 
temperature T. The field therefore presents random thermal fluctuations 
around its vanishing expectation value; they contribute to the energy Q and 
to the entropy. We shall study this problem by evaluating the free energy 
of the field, wherefrom all thermodynamic equilibrium properties follow; in 
particular its variations with the shape provide the constraints on the walls at 
fixed temperature. This free energy includes the zero-point energy, to which 
it reduces at T = 0. We shall thus treat simultaneously the Casimir effect 
proper, associated with the zero-point energy, that is, to virtual photons, and 
the radiation pressure effects for a black-body with arbitrary shape, which are 
associated with real photons in equilibrium with the walls that act as a 
thermal bath. 

Since general relativity is of no relevance in the present problem, energy 
is defined within an additive constant. We are interested only in its varia- 
tions and can thus get rid of divergences in the theory by substracting some 
constant that will tend to infinity. 

We wish to define the Casimir effect as a property belonging only to the 
field in vacuum, independently of the nature of the material in the walls. In 
general the zero-point energy of a field depends on the matter to which this 
field is coupled. We shall use the term "Casimir effect", in contrast to some 
authors, only when this energy can be defined separately. For electromag- 
netic fields within real materials, the interaction between the field and the 
charges does not allow us in general to separate out the energy Q of the 
field alone; moreover the presence of a material affects the field even out- 
side it. However, a complete decoupling is achieved in the limit of perfectly 
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conducting boundaries, which can be approached experimentally by use of 
superconductors. Both the electric and the magnetic fields E and B vanish 
inside them. Outside them, they simply impose the boundary conditions 

Et = , B„ = (3) 

on the tangential (t) and normal (n) components of the field. The presence 
of material bodies then has only a mere geometric effect. The part — J j.A 
of the energy, which involves the currents and the vector potential and which 
determines the matter-radiation coupling in the equations of motion, can 
be assigned to the matter and left aside while is assigned to the vacuum. 
Anyhow, for perfect conductors, in a gauge where E = —dA / dt, this coupling 
energy vanishes on average since A is perpendicular to the surface current j . 

In this idealized model, the field and the matter of the boundaries do 
not exchange any energy, even in time-dependent situations, although the 
coupling between the potentials and the charged particles relates E and B 
to the charge and current densities through Maxwell's equations of motion. 
Indeed, the rate of decrease of the energy Q in some empty region of space 
is the outgoing fiux of the Poynting vector /Uq ^E x B across the boundary of 
this region. The conditions imply that the Poynting vector is tangent to 
a perfectly conducting wall, and hence that no energy can fiow across such 
a wall. The establishment of thermal equilibrium in a vacuum is ensured 
only by the fact that real conductors are never perfect; this allows energy 
transfers between field and matter. Similarly to the model of an ideal gas, 
the present model is too crude to describe the establishment of equilibrium, 
but it is adequate for equilibrium properties. 

1.3 Synopsis 

We first recall Casimir's calculation for parallel plates, as well as its general- 
ization to arbitrary temperatures. Through the further use of the so-called 
Derjaguin approximation, this allows us to briefiy describe a recent experi- 
ment, and the irrelevance of temperature effects for the latter, allowing the 
conclusion that macroscopic quantum vacuum fluctuations are indeed ob- 
served at room temperature! 

In the further study of arbitrary conductor geometries, we shall exhibit 
general aspects of the Casimir effect, for walls with arbitrary shape and 
at arbitrary temperatures. We wish to answer a few theoretical questions. 
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Figure 1: Conducting plates configuration. 



Can the energy associated with the quantized electromagnetic field in the 
empty regions bounded by conducting walls be defined independently of the 
properties of these walls? How does it vary with the temperature of the 
photon gas that constitutes the field? Does the existence of the Casimir 
effect depend on the specific features of electromagnetism? 

We shall rely on two detailed articles jT3 CHj which deal with the above 
questions. This will allow us to leave aside the technicalities and to focus 
on the various ideas. The main result is embedded in section 15.51 For a 
bibliography on the Casimir effect we refer the reader to articles published 
in H , and to recent monographs [71[Ill2niinil2ai21l21l2Sl- 

2 Casimir's calculation (1948) 

Consider two identical plates, parallel to the plane yOz, with a large area 
A = Ly X Lz, and separated by a distance L from each other along the 
orthogonal direction Ox. 

The ideal plates are perfectly conducting, and the stationary modes are 
described by wave vectors {k^:, k), where k = {ky, k^) is parallel to the plates 
; in the perpendicular direction, the boundary conditions give the discrete 
series: k^ = — , where n G N, while periodic boundary conditions along the 
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plates give {ky,kz) = I — — , — — I , where ny,nz G Z. This eigenmode, 

\ Ly Lz J 

denoted by (n, k) , oscillates with a "frequency" uj 



^ = ^„(k) = cy-^ + k2. (4) 

Each eigenmode occurs with two possible polarizations (except for n — Qi). 

At zero temperature, the cavity electromagnetic energy is the sum £q of 
the zero-point energies of the eigenmodes 

£q= ^ £o['^n(k)] , £oM = (5) 

modes (n,k) 

The terms of this series are unbounded, and the series diverges. However, 
physics tells us how to cure this problem: at very high frequencies, the ma- 
terial ceases to conduct and becomes a dielectric, ultimately transparent to 
radiation. Then the boundary conditions no longer apply, and very high 
eigenmodes no longer contribute to the resulting force. This leads to intro- 
duce the mathematically convenient regularization 

S.^ E eoK(k)]xf^), (6) 

modes (n,k) 

where the cut-off function xi^/^c), is such that x(0) = 1, and is regular 
at the origin. It vanishes, along with all its derivatives, for u/oOc +oo, 
sufficiently fast so that the sum converges. The cut-off frequency ujc appears 
in X for dimensional reasons ; it depends on the microscopic characteristics 
of the material. (The perfect conductor limit corresponds to Uc +oo, so 
that xi^/^c) 1 for any finite u.) 

For large values of Ly, L^, one can replace the sums over parallel wave 
vectors by integrals. The distance L remains finite and associated with dis- 
crete modes. On thus gets 



modes (n,k) n=Q 



R2 



*As usual, the asymptotic result does not depend on these boundary conditions. 
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where the prime means that the n = mode has weight 1/2. Let us define 
e{u!) = Eq^uj) X i^/^c) ■ The energy ^ can then be written from ((7j) 



- oo „ 



Owing to (0} on has, for n fixed, the form uduj = c^kdk, where k = |k|. By 
simple integration over parallel wave vectors k: 

l> P+OO 

/ d^ke: [ti;„(k)] = 27rc^^ / ujdujE{uj), ujn{0) = ncn/L, (8) 



whence: 



n=0 ,j 
TTcn/L 



The associated force, Xq, can be readily derived, 

n=0 ^ ' 

The equivalent X,^ of Xq in the large L limit is given, as any continuum 
limit, by substituting the integral over n for the "primed" sum over n: 



oo 

r2i 



^^ = -^^ i^ng{n). (11) 

To obtain the zero-temperature force, associated with the vacuum energy 
£0(1^) = and acting on the plate, one must also take into account the 
opposite force exerted by the (infinite) electromagnetic vacuum outside of 
the capacitor. This force is just the opposite of (fTT|) . whence the resulting 
force 



2L^ 



00 00 
E '9{n) - Jdng{n) 

n=0 



To evaluate the difference between a series and the associated integral, we 
use the Euler-Maclaurin formula: 

00 +00 1 1 

E 'g{n) - J dn g{n) = --^'(0) + -g"'{0) + O {g^'\0)) , (12) 

n=0 J-Z 0! 



which involves all g^s derivatives of odd order, taken at the origin, and which 
is valid for a function g vansishing at infinity, as well as all its derivatives. 
By calculating the successive derivatives, one finds here 

^'(0) = 0, g"'{0) = 6x(0) = 6, ^7^(0) = O ,p>3. 

One therefore finds the finite value^ 

oo oo 1 

^'g{n)-jdng{n) = - + 0{u;,~'). 

71=0 

The zero-temperature resulting force thus possesses a universal limit for per- 
fect conductors, i.e, when +oo. The limit pressure, found by H. B. G. 
Casimir in 1948, is: 

= (13) 

The Casimir force is attractive, and one finds the analytic form which was 
anticipated in terms of c and of the length L. Only the numerical coefficient 
remained to be found: — 7r^/240, and the remarkable fact is that it is non- 
vanishing and universal, i.e., independent of the microscopic nature of the 
perfect conductors. To the resulting Casimir force is associated a subtracted 
zero-point energy Eq such that 

~A^'-'rmT^' ^^^^ 



3 Electromagnetic free energy of a planar ca- 
pacitor 

3.1 Vacuum and thermal parts 

Let us briefly consider the effect of temperature, in order to compare it to 
the zero-point effect^. At a given temperature T, photons will fill in the cav- 
ity, and they obey the "black-body" statistics. To each classical eigenmode 

^The presence of the factor , varying rapidly with n, yields the non- vanishing value 
(?"'(0) = 3!. In the absence of such a term, the Euler-Maclaurin formula would start with 
(7'(0) = O (wc""'^), and the difference between the sum and the integral would vanish in 
the perfect conductor limit. 

■tXhe first calculations are due to Fierz [5^ and Mehra j27| . 
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(n, k) is associated the Hamiltonian of a quantum harmonic oscillator, with 
frequency v = uj ^ lo = C(j„(k), given by (j3)). The eigenvalues of this Hamil- 
tonian are then e^^ = huj{N + 1/2), where > is the number of photons 
in the mode. 

The free energy of an eigenmode u at temperature T, with /5 = l/ksT, 
and ks Boltzmann's constant, has the form 

1 °° 11 
/(^) = _iin e-'^'^^^+s) = -nw + Ain (l - e"'''^) . (15) 

^ 7v=o 2/3 



We recast it as 
with 



fico) = eoico) + hico), (16) 

soiu) = hw, h{u) = iPhw) , (17) 

where /t(i^) is the thermal part of the mode free energy, and where 

(fi{x) = In (l - e-"") , ifi{x) < 0. 

By definition, the zero-temperature limit of /t(^^) vanishes. 

The purely thermal part of the electromagnetic free energy between the 
plates is then defined as follows: 

^T= /TK(k)]- (18) 

(n,k) modes 

In contrast with the eigenmode sum associated with the vacuum, the sum 
(fTH|) associated with the thermal radiation is convergent. By using ((7j) it can 
be written as 

W) = '2jf^,Y.' d^k [/3^„(k)] . (19) 
In (I19|l the integration over parallel vectors k gives, as in (jH)): 

!> i>+oo 

/ d^k<^[(]nuJr^{k)] = 2nc-^ udu <^{(3nuj), (20) 
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with ujn{0) = Tcnc/L. In terms of the dimensionless variable x = phu, and 
of 

+00 

■?/'(«) = J dx x(p{x), (21) 

u 

the free energy is obtained simple series: 

A I °° 

^^^^^ = 2^(^E>(^«)' « = ^^^^/^- (22) 
3.2 Continuous limit 

By comparing the energies hu of photons belonging to two consecutive eigen- 
modes, on can estimate the domain of temperatures T or separations L for 
which the discrete character of the eigenmodes disappears. For vanishing 
parallel wave vectors: hAu = h[LJn+i{0) — ujn{0)] = hnc/L = aksT, and 
the eigenmodes appear as a continuum for a < 1. At ordinary temperature, 
T ~ 300 K, this gives L > 24 /xm, and only for shorter distances will the 
discrete character of the eigenmodes be detectable. 

Let us introduce the large L limit, jFr^(L), of the free energy J-'t{L) (j22)). 
In the limit a = Pnhc/ L -C 1, the series in JF^ converges towards the integral^ 

+00 



The numerical coefficient is 

+00 +00 

du -ipiu) = dx x'^Lp{x) = —2^(4) = , 

J 45 


in terms of the Riemann ( function. The continuum free energy 



A 1 2^,., AL 1 



TT 



2 



§The presence of the "prime" notation in the sum over n and of the factor i for the 
n — mode are irrelevant in the large L limit. They imply a finite difference between the 
smns, which are here evaluated at order 0{L) as integrals. 
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is precisely the black-body free energy in a large volume Q = A x L. 
It is convenient to introduce the thermal subtracted free energy 



J-T — J-T — -^T > 

which reads, owing to and 



A 



5^V(«n) + -C(4) 



(24) 



(25) 



3.3 Thermal forces 

The force on, e.g., the right hand plate, is calculated as 

dJ^T{L) 



Xt{L) = 

From the free energy (j22j) one finds 



dL 



(26) 



A 



27x131 {phcy 



(27) 



n=0 



This force is perpendicular to the plate and positive {(p{x) < 0), hence re- 
pulsive. It is the black-body pressure in a finite geometry. In the L oo 
limit, one recovers the infinite volume black-body radiation pressure 



A 1 



(3 {(3hcfAh' 



(28) 



In fact, a given plate at temperature T is in equilibrium with the thermal 
radiation existing on both sides. It thus will also be subjected to the pressure 
force exerted by the exterior photons, which is just the infinite volume black- 
body pressure force, — X-^, that we just evaluated. The resulting thermal 
force is therefore: 



x^ 



d (J^T - ) _ d^T 



dL 



dL 



(29) 
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3.4 Short distance or low temperature expansion 



The series ()27|) giving the interior thermal force Xt yields a natural low- 
temperature or short- distance expansion for a = jSirhc/L ^ 1. Indeed the 
n = term vanishes, and for n > 1, ip{na) ~ — e~"". The n = 1 mode thus 
contributes a leading exponentially small repulsive force: 

Xt = [e-^ + O (e-2")] , a » 1. 

In the thermal resulting force ()29p Xt, the leading term will thus be the 
black-body exterior term: 

3.5 Comparison of zero-point and thermal effects 

Adding (fT3|l and (|30|l . one gets the total pressure at low temperature or short 
distance: 

= h + = --kM ~ ^^\wcf ^ w'''^ + ■ ■ ■ ■ ^^^^ 

Thus the resulting force is dominated by the Casimir and black-body forces, 

both attractive. For L = Ifim, indeed one gets a ~ 24. Hence the first 

internal mode contributes to the thermal force only with a relative factor e~^^! 

Xx X°° 1 /27r\^ 
The ratio 7 = is thus, owing to (IHTl) : 7 ^ ^ = - I — I . For L = 

Xo _ Xo 3 V a / 

500 nm, a = 48, one gets 7 = 0.98 x 10~^. Thus even at room termperature, 

the vacuum fluctuations largely dominate the black-body effects. This is 

equivalent to a zero-temperature situation, and sensitive experiments will be 

able to detect the quantum vacuum effects. 



3.6 Total free energy 

Let us finally introduce the total free energy F associated with the vacuum 
energy and the thermal free energy^ 

F = Sq + ^t = So + J^t- T^. (32) 

^In this total free energy the contribution of the exterior volume has been subtracted 
out. Its definition coincides with that of the regularized, then renormalized free energy 
given in sections l5 . 21 and [5?5l below . One has in particular: F(0) = F{T) — F(0) = J-^. 
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The total force acting on a plate is therefore 



dF 



(33) 



Using (fT^ and (|^. we have 



A 



L3 



720 



Q{a) = - 



X;V(«^) + -C(4) 



n=0 



(34) 
(35) 



The above series yields the natural low-temperature expansion of the thermal 
function Q is for a ^ 1. Owing to definition 



V'(O) = -C(3), ^(a) = -(a + 1) [e-" + O (e"^-)] , a > 1, 



(36) 



whence 



L3 



1 



1 



720 



-^C(3) + ^C(4) - (a + 1) [e- + O (e"^'^)] 



(37) 

for a = p-Jihc/L ^ 1. 

The high-temperature expansion can be obtained from (jH^ and the Pois- 
son formula. One can also use a remarkable duality formula between low and 
high temperatures H ^S] 

a'^g{a) = a'^Gia'), aa' = {2tiY, (38) 
which yields for the thermal free energy: 

:Ft{L, a) = (27r/a)^ {L, (27r) Va) • (39) 
From (|37|) we deduce the high-temperature limit of the total free energy: 



+ O (/^-^e-^'^'/") , a = (3nnc/L « 1. 



(40) 



"This duality is valid only for parallel plates, see section W>A\ below. 



14 



o 



L 

O' 



Figure 2: Experimental set-up of the conducting surfaces. 

4 Force between a sphere and a plane 

4.1 Experimental situation 

An actual experimental set-up is described in Figure |21 A metallized sphere 
of radius R is placed in front of a conducting plate, at a distance 00' = L. 
In the experiments, like the one performed in 1998 [H], a polystyrene sphere 
is attached to the arm of an atomic force microscope, and placed in front of 
a polished planar surface. Their surfaces are coated with an aluminium layer 
a few hundred nanometers thick. To prevent corrosion, they are additionally 
coated with a very thin alloy layer, which is transparent to the radiation. The 
overall radius of the sphere is R = 98.0 ± 0.25 /im. The range of distances 
is 120 nm < L < 500 nm. Measures are performed at room temperature. 
Measurements with parallel plates, as in the original Casimir's calculation, 
are also performed jTT], but are more difficult, due to the necessity to properly 
align the plates. This experimental problem is avoided in the sphere-plate 
geometry. However, the theoretical force is not explicitly known in this case, 
eventhough it exists in a closed form (see sections 1^751 and IB. 5. 31 below). 

4.2 Derjaguin approximation 

An approximation method, due to Derjaguin (1934) j^H], allows in the L <ti R 
limit, the calculation of the sphere-plane interaction in terms of the purely 
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planar interaction. One replaces each elementary slice of the sphere cut 
parallel to the plane, by its orthogonal projection towards the plane. The 
resulting force then is 



A 



Xix) dx, 



(41) 



where X{x) is the force between two planes at distance x. From ()41j] 
and definition of X we immediately find 



X^P'^(L) = 2nR-F{L), 



(42) 



where we used the fact that the total free energy F of two plates vanishes at 
infinity, as shown by the large-distance or high-temperature equivalent ()4U|). 
The Casimir zero-point force on the sphere is thus: 



xf^{V) 



vr^ he 
360 Z^' 



(43) 



where is the Casimir energy ()14|1 of two plates. For R = 98.0 ± 0.25 /im, 
and for L = 200 nm, one finds for example: Xq 



sph 



-33.4 ±0.09 X 10"^^ N. 



Such a force of the order of tens of pico-newtons is perfectly measurable, and 
comparable to forces implied in biological systems, e.g., in micromanipula- 
tions of single DNA macromolecules. In the experimental range 120 nm < 
L < 500 nm, the minimal value of (L = 500 nm) = 48. Hence 

we can use the short-distance or low-temperature estimate ()37|) of the free 
energy F: 



""'-R^ < 1 



360 L3 



720 



2a3 



C(3) - ^C(4) + O (a-'e--) 



■ (44) 



The thermal correction in is governed by the C(3) term, which gives in 
the experimental range (L < 500 nm) a relative thermal correction 4 x 10~^. 
It is attractive and adds to the leading zero-point term (jlHj) . Formula (jiijl 
is used by experimenters jH]. 

The simple zero-point formula ()43|) is shown in figure El in comparison 
with experimental results. The match is good, eventhough the experimental 
curve goes further above (j43p at low values of L, with a relative shift in the 
tens of percents. This cannot be accounted for by temperature corrections 
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Figure 3: Comparison of the zero-point force (j43p (dashed) with experimental 
results; the full curve takes into account finite conductivity and rugosity 
corrections (ref. jB]). 

of order 10^'^, which furthermore lower the theoretical predictions. Good 
agreement (~ 1%) is found when finite conductivity and rugosity corrections 
are taken into account [HI El ISO] ■ Let us also remark that the analysis of 
geometrical corrections to the Derjaguin approximation (PT|) . following the 
general formalism described in the next sections, would be useful. 

Before addressing the Casimir effect for arbitrary geometries, let us con- 
clude that one does observe, even at room temperature, macroscopic electro- 
magnetic forces, generated by vacuum fluctuations, and proportional to he, 
and in the absence of charge and photons in the cavity! ** Notice that sev- 
eral experiments also give clear evidence of retardation effects in atom-wall 
interactions, in agreement with the Casimir- Polder prediction pT] . 

**It is interesting to note that there is an infinity of longitudinal soft photons (n = 0, k 7^ 
0) present in the cavity, but they do not contribute to the force; in contrast, photons with 
n > 1 are essentially absent from the cavity at such short separations L, which allows a 
direct observation of the vacuum energy, even at room temperature. 
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5 Arbitrary conductor geometries 



5.1 Eigenmodes 

Let us first consider a single connected region v limited by perfectly conduct- 
ing boundaries. The electromagnetic field in v can be analyzed in terms of 
the eigenmodes m, obtained by solving Maxwell's equations 

curl = iuJmBrn , curl B„ = -iu;Em , 
div E^ = , div = 0, 

with the boundary conditions We shall keep aside the electrostatic and 
magnetostatic solutions with zero frequency, which do not contribute to the 
Casimir effect. Each mode m behaves as a harmonic oscillator with frequency 
= to'm/STr = where has the dimension of an inverse wavelength. 

Its associated energy Q may take the quantized values ^(gm, ^) = ^cgm(-^+ 
|). At finite temperature it yields a contribution /(gm) ffTH|l to the free 
energy, which can also be written as (hereafter we use temperature units 
where ks = 1) 

/(g„) = rin[2 sinh(ncg^/2r)] . (46) 

At zero temperature the corresponding contribution to the energy of vacuum 
is \hcqm, the limit of as T — > 0. For the high frequency modes such that 
hv ^ T, the free energy /(g) is dominated by this zero-point energy. For 
the low frequency modes such that hv <^ T, 

/(g)~-Tln(r/;icg) (47) 

is dominated by the classical — TlnT behavior. 

The spectrum of eigenfrequencies, or equivalently of eigenwavenumbers 
in the considered region v is characterized by the density of modes 

p^^\<i) = Y.^{q-& , (48) 

m 

and the free energy of this region is formally equal to 

/>oo 

F(^)= / cigp('')(g)/(g) . (49) 



A first difficulty arises if the domain v is infinite, since the spectrum is 
then continuous. We therefore imagine that the full system is enclosed in 
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a large box S, with volume V, which will eventually tend to infinity. By 
assuming this outermost boundary to be perfectly conducting, we not only 
discretize the spectrum m in the open regions v, but also confine the field 
and ensure that no energy is radiated outwards. 

A second difficulty is associated with the fact that the spectra are not 
bounded. Indeed, for large q, the distribution (j48|) has the asymptotic expan- 
sion 



+ 



O (^^^ + 0{q'/^ X osc) . 



The dominant term is proportional to the volume v of the considered region; 
it is the only one which contributes to the black-body radiation in the ther- 
modynamic limit. The second one is a curvature term; it is the integral over 
the boundaries of v of the average curvature R^^ = ^{Ri^ +-R^^), where Ri 
and R2 are the two main curvature radii at the point a, oriented towards the 
interior of v. It is supplemented, in case the boundary has not a finite curva- 
ture everywhere and includes wedges with a dihedral angle 6 at the point s 
of the edge, by the next, wedge term, integrated along the wedge; we have for 
instance = 7r/2 if f is the interior of a cube, = 37r/2 if it is the exterior. 
Even if it is smoothed, the distribution (fSUj) finally includes oscillatory terms 
with an amplitude which increases with g^"!". Since f{q) ~ ^hcq for q ^ 00, 
all the terms exhibited in (jSUj) lead to divergences in the free energy 



5.2 Regularization of the free energy 

We encounter here the simplest example of the divergences that plague quan- 
tum field theory. We deal with them by using the standard technique. We 
first regularize the divergent formulae by means of cut-offs. We then deduce, 
from the resulting finite expressions, quantities that are physically observ- 
able at least theoretically. We finally renormalize these quantities by letting 
therein the cut-off parameters go to infinity. The theory is renormalizable if 
we get a finite limit for the physical quantities. 

We have already regularized the ^Hnfrared^ divergence associated with 
the infinite size of vacuum by introducing the box S. We shall deal with the 

''"''The corresponding semiclassical expansion of the density of modes over periodic orbits 
|17j has been recently experimentaUy tested in a superconducting microvawe cavity [IV2\ 
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^''ultraviolet divergence associated with the high frequencies eg ^ oo in the 
integral by introducing, as in section |21 above, a cut-off factor xiq) close 
to 1 for g -C Q and decreasing sufficiently fast for g 3> Q so as to restore 
convergence of the integral. Our final goal is to construct a renormalized 
free energy F, finite in the limit as E ^ oo, Q oo. If this is feasible, 
it will mean that the ideal model of the electromagnetic field outside a set 
of perfectly conducting boundaries is renormalizable. In other words the 
Casimir effect exists as a property of the field proper, conditioned by the sole 
geometry of the boundaries S. 

Let us first see how one can get rid of the most severe divergence, asso- 
ciated with the first term of (j5(jp which after regularization with x(g) yields 
F'^'"^ oc vhcQ^. Consider, for instance, the Casimir force between two spheres 
with volumes vi and V2- There is here a single empty domain vo, which lies 
outside the spheres and inside S. Its volume is fo = ^ — "J^i — "^^2- The only 
universal and natural way to cancel the corresponding divergence consists in 
replacing the solid conductors by hollow thin conducting shells, and in taking 
as a reference the free energy of the empty enclosure S, which is itself di- 
vergent. We now have three empty domains, Vq, Vi, V2, the volumes of which 
sum up to V. The most divergent term thus disappears if we substract the 
free energy of the empty space within E from the total free energy in the 
presence of the two spherical shells, as 

p{vo) ^ _^ p{v2) _ ^ ^^Y^ 

More generally and more precisely (see Fig. HJ, we denote as S the set 
of two-dimensional surfaces which bound the considered conductors. They 
partition the whole space (within the enclosure S) into a set of connected 
regions f , some of which coincide with the actual vacuum (as vq above), the 
other ones with the interiors of the conductors (as vi and f 2 above) . We then 
define the regularized free energy associated with the whole space partitioned 
by S as 



Freg = / dq 



f{q)x{q) = dq 5p{q)f{q)x{q) ■ 



^p('^)(g)-p(^)(g) jy^j^y^j-i ^..^w;.w;aw. 

(52) 

This expression is finite owing to the cut-offs S and x(Q')- If; indicated 
above, it has a finite limit F as S — 00 and x(Q') independently 
of the shapes of S and x{(l)y the Casimir effect will appear as a universal 
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Figure 4: Set S of surfaces bounding the interiors Vi and V2 of conductors, 
and defining a third connected (vacuum) domain vq within enclosure S. 

property characterized by the free energy F for the boundaries S, which will 
depend only on the geometry of S and on the temperature. This will provide 
us with a generating function for all mechanical and thermal properties in 
thermodynamic equilibrium, in two idealized circumstances. 

On the one hand, the expression can be directly interpreted as the 
change in the free energy of the vacuum when a system S of closed, extremely 
thin conducting foils is introduced. The variations of F under deformations 
of such foils determine the constraints induced on them by virtual (for T = 0) 
or real photons (for T 7^ 0). For instance, for a single sphere separating two 
regions f 1 and V2, the dependence of on the radius determines the pres- 
sure exerted on the skin of this hollow sphere by the internal and the external 
field. We shall also encounter below constraints which tend to corrugate such 
thin sheets, by studying how F changes under periodic deformations. 

On the other hand, the expression is also suited for the study of forces 
between bulky indeformable conductors, as we now show. Let us return to 
the above example of two spheres. After regularization the force between 
them is associated with the variation of F^'"°'' when they are shifted apart. 
The regularized quantities F^^^^^ , F^^^^ which enter fj51|) are not physically 



21 



relevant to the present problem where we deal with bulky rather than empty 
spheres, but they do not depend on the distance between these two spheres. 
Actually, a perfectly conducting skin of a sphere behaves as a perfect screen 
and the electromagnetic fields, inside and outside, are independent. Thus the 
force evaluated from ()51|) is the same as that evaluated from F^'"°\ and it is 
preferable to use because the divergences are expected to be eliminated 
by this combination. More generally, whenever solid conducting bodies can 
be displaced but not deformed, we can derive the forces between them from 
the Casimir free energy ()52j) for which the interior of each body is replaced 
by vacuum. This trick will allow us to renormalize F. (However, for thermal 
properties, one should leave aside the contributions of real photons within 
the conductors.) 

As mentioned above, the cut-off factor which regularizes the integral 
fl52|) for large q has a physical meaning. At high frequency, real conductors 
are never perfect. Electromagnetic waves can penetrate them, and go freely 
across them if they are thin. The objects S become transparent and the 
modes within S tend to be the same, whether S is present or absent. Thus, 
for imperfect conductors, the factor Sp{q) in would decrease for large q. 
In our model we simulate imperfect conduction or transparency at high fre- 
quency by evaluating Sp{q) for perfectly conducting sheets S and multiplying 
by 

We now proceed and study the behaviour of (|52j) when the boundary E 
is pushed away to infinity and when the conducting sheets S tend to become 
perfect with — >• 1. 

5.3 Fields in the presence of perfect conductors 

Our strategy will rely on the following ideas. 

(i) We replace the solution of eqs. ©, (03) , which define the modes in 
each region v, by the determination of the Green functions associated 
with these partial differential equations and boundary conditions. Such 
a Green function contains in a synthetic way the whole information on 
the modes. It is a function of a complex variable k, analytic in the 
upper half-plane. 

(ii) We express the distribution of modes p^^\q) in terms of the boundary 
value for A; — 5> g + iO of the Green functions. 
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(iii) This will allow us to regard (|52p as an integral in the complex plane k 
along the half-line k = q + iO, q > 0, and to deform this contour towards 
the pure imaginary axis k = iy, y > where the Green functions are 
more regular than along the real axis (they have an infinity of poles at 
k = ±qm)- 

(iv) We determine the Green functions by means of Neumann's method, 
which expresses them as solutions of two-dimensional integral equations 
over the boundaries S and S. The regularized free energy will thereby 
be expressed through the kernel of these integral equations in terms of 
the geometry of S and S. 

(v) Along the new integration contour k = iy, we can solve these integral 
equations by iteration. The resulting series are convergent, and can be 
interpreted physically as describing multiple scattering of an electro- 
magnetic wave on the walls, involving successively induced currents. 

(vi) Convergence of the multiple scattering expansion allows us finally to 
control the limit S — oo, % ^ 1 and to find an explicit expression for 
the limit F of (j^ . 

We shall content ourselves here with a sketch of this programme. Detailed 
proofs can be found in pTl llSj. Given the symmetry between the fields E and 
B in eqs. (II3|) . it is convenient to introduce two Green functions, a magnetic 
one r(r, r') and an electric one S(r, r'), which are tensors with two indices at 
r and r'. The first one represents the magnetic field created at the point r, 
in the presence of the conducting boundaries S and S, by a magnetic dipole 
lying at r' and oscillating as e~'^'^* at the complex frequency ck/2TT. The 
current density associated with this source is (within the factor /io) 

jo(r, r') = curl^ [5^{r - r')l] (53) 

where 1 is the unit tensor. In each region v,r{r,r'; k) is the solution of the 
partial differential equation 

(V^ + fc^)r = -curl jo , div r = 0, (54) 

obtained by eliminating E from the eqs. ()45|) . with the boundary conditions 

Tn = , (curl r)t = 0. (55) 
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The Green function r(r, r') is the magnetic field generated at the point r by 
the source jo at r' and the currents j(a,r') that it induces at the points a 
of the conducting surfaces S*, S. We shall denote as Ua the normal vector 
at a, oriented towards the region v where r and r' lie. Both jo and j are 
tensors with two indices; the first one refers to the direction of the current, 
the second one to the orientation of the dipole at r'. Using the formalism 
of retarded potentials, we can express the magnetic field created by each 
elementary current j as Mj, where M is the kernel 

M(r,r') = curl, [Go(|r - r'|)l] . (56) 

A product like Mj stands for integration over space and summation over a 
tensor index. The scalar Green function 

"■.(O = (57) 

is the solution of the equation (V^ + k'^)GQ{r) = —6^{r) that vanishes at 
infinity for Imk > 0. The magnetic Green function F is thus expressed as 

T{r,r') = [ £r"M{r,r") }o{r",r') + [ d^a M(r, a) j(a, r') ; (58) 

its first term To is the field produced by the dipole in the infinite space. 
The as yet unkown currents j satisfy the integral equation on the boundary 

j=ji + ^j , ji(«,r) = J d\K{a,r)io{r,r'), (59) 

where the kernel K between two points a and /5 of the boundary is the tensor 

K{a, (3) = 2n„ A curU[Go(|a - (60) 

In eq.()59|) the product Kj stands for summation on the tensor index and 
integration over a on the boundary. The proof of (j59j) . (j6(J|) relies on an 
extension of Neumann's method, based upon the discontinuity of the surface 
integral in (j^Hjl when r crosses the boundary p[B1. 

Altogether the solution of ()59|) determines j, and F follows from ()58|) . 
taking into account the various definitions (j^Hj) . (j^ . ^7\i . (pljl . The solution 
of the partial differential equation (j^^ with the boundary conditions ()55p 
thus amounts to the solution of the integral equation (j59|) on the boundary. 



24 



The Green function F is a generating function for the modes m defined by 
©, (03) in the connected region v where the source hes. Indeed, in terms 
of the complex variable fc, its poles are the real points k = ±qm = ^^^m/c 
and the corresponding residues are given by 

r(r, r'; k) = J2 ^^B„(r) ® B„(r') , (61) 

m 

where the magnetic field Bm{r) for the mode m is the real solution of (jHl), 
()45|1 ■ normalized according to J (i^rB^(r) = 1. The spectral density pHjl is 
thus related to T through 

p('^) (q) = — [ d^r tr Im T{r, r;q + iO) , (62) 

where the trace tr refers to the tensor indices. 

We can likewise introduce an electric Green function H by interchanging 
magnetic and electric fields, which amounts to interchanging the boundary 
conditions (j55j) . It is the electric field created at the point r by a source with 
current density curljo/icfc in the presence of the boundaries S, S. Taking (j53|) 
into account, we see that, except at r = r', this source produces the same 
electric field as an electric dipole with current density —ik[6^{r — r')l]/c. 

Provided r and r' lie in the same region v, the electric Green function 
S(r, r'; k) can be represented by an expression analogous to (jKHjl . (jSHj), within 
the mere change in sign o/ji. (However, whereas T is expressed by means 
of ()58|1 in terms of a true electric current density, the similar expression for 
H involves a fictitious current, without physical meaning. Moreover, if we 
take r and r' on different sides of a boundary S, the expression (j^Hj) for T 
vanishes as it should, whereas the similar representation for 3 provides an 
unphysical non-zero value. Such a behaviour is currently found in the books 
of mathematics that deal with Neumann's method; in fact, paradoxically, 
single-layer potentials are used there to represent Green functions when they 
vanish on the boundary, double-layer potentials when their normal derivative 
vanishes, whereas it is natural in electrostatics to make the converse choice. It 
is the application of Neumann's method in two different ways which allowed 
us to find the above simple relation between the representations of F and H.) 

Since H has the same spectral representation ()61|) as F within the re- 
placement of Bm by E„i, we can equivalently express the distribution of 
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eigenmodes (jU^ as 

[q) = — [ d^r tr Im[r + H(r, r; g + iO)] , (63) 

and simplications will appear owing to this combination. 

In fact, the iteration of the integral equation (j59|l provides 

j=ji + A'ji + ir2ji + ... , (64) 

a series which exhibits the surface current j as the sum of successively in- 
duced currents: ji is according to a current directly induced by the 
dipolar source on the conducting boundaries; it induces in turn through the 
propagator K a secondary current K]i, and so on. The expression (jUUj) of K 
and the behaviour of the free Green function (j37jl show that this propagator 
decreases exponentially at large distances (while oscillating) for Imk > 0. 
Moreover, K{a,f3) vanishes when the point f3 lies in the plane tangent at 
a to the boundary. Thus a current circulating on a plane boundary does 
not induce through K any secondary current on the same plane. At short 
distances K vanishes for a smooth boundary and is proportional to its cur- 
vature; this ensures the convergence of the integrals in Kji, K^ji, etc. By 
relying on these properties, we can show that the expansion ()64|) is conver- 
gent at least in the region Imk > |Re/c|. (For k = describing static fields, 
the convergence depends on the topology of the boundaries.) The general 
theory of the Casimir effect will make use of this convergence. 
The series for F which results from (j^Hj) and reads 

r = M--^ 30 = M 30 + MK 30 + MK^ 30 + MK^ jo + ■■■ , (65) 
i — A 

where M defined by describes the propagation of a wave issued from 
a unit element of current, and where K defined by (jUUI) describes a similar 
propagation followed by the creation of an induced current. We can thus 
interpret ()65j) as a multiple scattering expansion: the wave issued from the 
source at r' may reach directly r (first term); it may propagate from r' to 
a point a of the boundary where it is scattered to reach r (second term); 
it may scatter successively twice on the boundary before reaching r (third 
term); and so on. The grazing scatterings vanish, so that the series fl65p 
reduces to its first two terms for a single plane boundary, the second one 
describing the reflected wave. 
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The electric and magnetic Green functions are related to each other 
through 

H(r, r'; k) = k~'^ curl^ curl^' [r(r, r'; k) — r(r, r'; 0)], (66) 
and conversely. However the alternate series 

H = M-^ jo = Mjo-Mir jo + Mir^jo , (67) 

which results from the integral equation over 5, S for S does not correspond 
term by term to (jMj) . 

This expansion has the same interpretation in terms of successive scatter- 
ings as (j65|) in the physical situation when r and r' lie in the same connected 
region v. (However, if r and r' lie in two different regions v separated by 
a boundary the series converges towards some non-zero value, and 
thus does not represent the actual electric field produced at r by the source 
at r' . On the contrary the various terms of the expansion interfere 
destructively in such a configuration.) 

Adding (jHKjl and (jHTj) as in (jMjl cancels all the odd terms of the expansion 
and yields a geometric series in K^. Only survive in the evaluation of the 
distribution of modes the terms describing an even number of scatterings. 
On the other hand, in the summation over the regions f , we shall have to 
transfer the two points r and r' (with r = r') from one side to the other 
of the surface 5*. The kernel K changes its sign in this operation, since its 
definition ()60p involves the normal vector Ua oriented in the direction of the 
domain v where the Green function is evaluated. After addition of T and 
H in the eigenmodes on both sides of 5* are evaluated with the same 
integrand. This remark will allow us to perform explicitly the summation 
over V and the integration over r in (j68|) below. 



5.4 Limiting process 

We are now in position to express the regularized free energy ()52|) in terms 
of the kernel K on the boundaries S, S. In order to switch the integration 
over q towards the complex plane fc, we introduce the generating function of 
the modes: 



6^{k) = - [ d^r tr lim 

2 J r'—tr 



(rM+H(''))-r(^)-H(^)(r,r';A;) 



(68) 
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As in (I51|) or (j52j) . we have summed over the various regions v bounded by 
S and enclosed in S, and substracted the contribution of the empty box; the 
integral is therefore carried over the whole interior of S. The Green functions 
are singular for r' —>■ r, but their singular part, which arises only from the 
first term To{r — r') = So(r — r') = Mjo of (jHB) and ((Hf ), is cancelled in (|nH|l 
by the subtraction of T^'^^ + H^^^. Thus no divergence appears in (j68|) . 

Like the T and H's, the function 5$(/c) has no other singularity than poles 
at the points k = ±qm of the real axis, with residues =F|'?m- Hence, according 
to (jSnj); we have 

Sp(q) = —Im 6<^(q + iO) (69) 
Tvq 

for g > 0. The expression (jU^ of the regularized free energy thus 

reads 

Freg = — Im dk ^ ^ ^ ^ In (^2 sinh— J xik) . (70) 

The subtraction of the real number (5$(iO), which is equal to 5$(0) when the 
box S is finite, ensures the convergence of the integral; it corresponds to the 
fact that static fields associated with A; = do not contribute to the Casimir 
effect. 

In order to take advantage of the convergence of the expansions and 
(jnZj) for Imk > \Rek\ in (jHEI), (HOI), we shall deform the integration contour 
in (fTOI) towards the imaginary axis. We also get rid of the spatial cutoff E. 
We therefore introduce the function 

^(y) = lim 5^{\y) . (71) 

When taking the limit S oo, we note that Go defined by (j37j) and hence 
M and K defined by ()56|) and ()60|) decrease exponentially with the distance 
/ as e~^'. The function $(i|/) is represented, through its definition (j68|) and 
the multiple scattering expansions (j65p . (j67|) . by closed paths that bounce an 
even number of times on or E. All paths which involve only scatterings on 
S are compensated for, owing to the subtraction in ()68p. The contributions 
of all the paths which involve at least one back and forth travel between 5* 
and S contain a factor e~^^^ where L is the minimum distance from S to 
E; they disappear when S is pushed away to infinity. The remaining paths 
which contribute to '^{y) involve only scatterings on 5, in even number; 
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before integration over r' = r of (j68p they contain for large r a factor e~^^'^, 
where L is the shortest distance between r and S", so that the integral over r 
is convergent. We shall perform it explicitly below. Note that the poles on 
the real axis of [5^{k) — 5$(i0)]//c with residues ±|, above which the integral 
(fTOj) runs, become dense in the limit S ^ oo and are replaced by a cut. 

It remains to get rid of the ultraviolet cutoff x{k) which eliminates the 
high frequencies. In order to deform the contour of (j7Up towards the half line 
k = iy,y > where the singularities of the logarithm lie, we take for x(^) a 
meromorphic function of the form 

The poles /Xj lie in the first quadrant and their residues satisfy 

^ Re ai = , -2 XI I^e(ai//i- ) = 1 . (73) 

i i 

The moduli have the same order of magnitude as a number Q that will 
tend to infinity, so that xil) is close to 1 as long as g -C Q and decreases as 

when q ^ Q. This behaviour ensures the convergence of (|52|) on account 
of the behaviour of the residual terms in ()50p . 

Deforming the contour of ()70|) towards the imaginary axis produces terms 
associated with the residues at the points k = fii. The properties of 5$(fc) 
then ensure that these residues vanish when oc Q — >■ oo. 

Indeed, for — > oo in the first quadrant, it turns out that the multiple 
scattering expansion (|^. (jU7|) is dominated by its lowest, two-scattering 
term MK'^jo, provided the surface S is twice differentiable. The result is 
found from a short-distance expansion as 

*'^'-3^//^"(«;k-l)^K7)' '''' 

where Ri and R2 are the two principal curvature radii and 1/R the average 
curvature at the point a. Finally, near the origin, \l/ behaves as 

^{y)^-n{l-Ay)+0{y^) , (75) 

where n is the genus of S, depending only on its topology (n = for a sphere, 
n = 1 for a torus), and where A > 0. 
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Figure 5: Sawteeth function g (I77j) drawn as a function of the dimensionless 
variable y/rj = yhc/27iT. 

5.5 The renormalized casimir free energy 

Altogether, using the above properties, we find the limit of (j52j) for S ^ oo 
and Q oo or — >• 1 as 

fir r°° dv 

F=- dy[^{y) - ^{oo)] + 2T / -^[^{y) - ^{+OMy) , (76) 

Jo Jo y 

where the temperature appears through the sawteeth function 

9{y) = ^- - + '^0{y -nri) , ri = . (77) 

^ n=l ^ 

with 6{x) = for x < 0, d{x) = 1 for a; > 0. A closed expression for the 
function \E'(?/), which encapsulates the effect of the geometry of the boundaries 
S on the modes of the field, is found by integrating on r in the whole 
space. We noted that its integrand is the same for all the regions v separated 
by S, so that this integral depends only on the first and last scattering 
points on S in the expansions ()65|). ()67p. The integration over r = r' of the 
corresponding product of M and Kjo can be expressed only in terms of the 
kernel K on S, and it simply yields ^y dK/dy. Hence we find 
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Figure 6: Closed paths (all with an even number of segments) contributing to 
the operator trace Trln(l — K'^) in the function ()78|) . Multiple scattering 
on a single shell contributes to the self-energy, while scattering between two 
surfaces contributes to their mutual energy. The enclosure S (dotted line) 
has been taken to oo. 

^(l/) = -|^Trln(l-ir2) , (78) 

where the trace Tr and the products stand for integration over S* of a variable 
a and summation on the tensor indices of K (see Fig. Expressed as 
function of y, the kernel K on the surface S defined by (j57jl . (j60|l reads 

-^-y\a~l3\ 

_\a-p\ 

It is real, decreases exponentially at large distances, and locally vanishes 
as the product of the distance and the curvature for |a — /5| — > 0. When S 
consists of several disconnected pieces, the normals Ua on each of them should 
be oriented compatibly; for instance, all of them should point towards the 
outermost region. Along the integration path of (fTBjl . the series obtained by 
expanding ln(l — K"^) in powers of converges. The values of "^{y) at both 
ends of the integration path are given by (f7i|) and (f7H|) . 



K{a, (3; y) = A curl^ 

ZTT 



(79) 
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The two terms of (|76p correspond to the two different phenomena that 
we are studying. The first one is the Casimir energy proper, associated with 
the variation of the zero-point energy of the electromagnetic modes that is 
induced by the introduction of the perfectly conducting shells. It is the 
product of he by a factor with dimension depending on the shape of S. 
The second one is the variation of the free energy of the black-body due to 
the effect of the boundaries on the gas of real photons. It has no ultraviolet 
divergence. The temperature T is that of the walls, which carry random 
currents in equilibrium with the quantized field. 

When there are several material bodies, the surface S involves several 
disconnected sheets. One can then classify the various terms of the expan- 
sion of (|78|) in powers of K according to the position on these sheets of the 
scattering points ■ ■ ■ (Fig. The terms for which all these points lie 
on the same connected sheet describe the free energy of each separate body; 
they do not contribute to the forces between indeformable bodies but de- 
termine internal constraints for thin conducting foils. Those for which some 
propagation K{a, (3) occurs between two points a and (3 situated on different 
sheets describe the interaction free energy. 

5.6 Conditions for the existence of the Casimir effect 

We have proved above the existence of a limit F for the free energy of the field 
only, for perfectly conducting, thin walls, which defines the Casimir effect. 
This has been made possible under two conditions, which we now discuss. 

On the one hand, the Casimir effect exists only owing to the electro- 
magnetic nature of the field. Our proof made use of the cancellation of the 
one-scattering terms iM/^jo of the expansions (j65|) and (j67|) . which elim- 
inated a divergent surface contribution. Let us show that the ultraviolet 
divergence of cannot be removed for a scalar field, so that the specific 
features of electromagnetism are essential to assign an energy to the field sep- 
arately. Actually, the high-frequency expansion of the density of eigenmodes 
that replaces (|5(J|) for a scalar field includes, after a volume term fg^/27r^, an 
area term equal to — sg/Svr for Dirichlet boundary conditions (cancellation 
of the field at the wall), or to +sq/8n for Neumann conditions (cancellation 
of the normal derivative), where s is the area of the boundary of the domain 
V. This term is associated with the occurence of single scattering in the ex- 
pansion analogous to (j65j) . In the evaluation of the regularized free energy 
(j52|) the contributions of the two sides of the walls add up. A contribution 
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^{hcs/Sir) dq q'^xi^) proportional to the area s of the walls S occurs, 
and its divergence when the ultraviolet cutoff Q tends to infinity is incur- 
able. A divergent Casimir force, tending either to stretch or to shrink the 
boundaries, would thus appear for a scalar field. The properties of matter 
interacting with the field could not be disregarded. 

The inexistence of a divergent area term for the electromagnetic field can 
be traced to the mixed boundary conditions Et = 0, (curlE)n = and to 
the constraint divE = 0. The first condition is of the Dirichlet type for the 
two tangential components of E; the constraint yields a Neumann boundary 
condition for the normal component. However the condition (cur IE) „ = 
relates the two tangential components to each other, so that altogether the 
electromagnetic field behaves as two scalar fields, one with Dirichlet, the 
other with Neumann boundary conditions. The fact that the area terms are 
opposite for these two boundary conditions entails their compensation. 

A second condition is also necessary for the renormalization of the total 
Casimir energy, namely the smoothness of the boundaries S. We have relied, 
in our elimination of the ultraviolet divergence, on the behaviour (j74j) of "^{y) 
for large y, which itself requires a finite curvature of 5*. To understand the 
origin of this condition, let us return to the expansion ()50p of the eigenmode 
density p^""^ for a given volume v. Its first term was cancelled by the sub- 
traction p*^^) in ()52j). The next curvature term of (jSUj) yields in F^^\ given 
by (jiUI) . a / qdq divergence. Fortunately, at each point a of the boundary 
5", the curvatures are opposite on the two sides of this boundary. Hence, 
in the summation over the domains v of (jH2I); the curvature terms cancel 
one another. Accordingly, the large y expansion ()74p of "^{y) begins with a 
second order curvature contribution. 

If, however, the surface S has a sharp fold, for instance if it is a hollow 
thin cube, this cancellation of the divergences from both sides 5* no longer 
occurs. The third, wedge term of (jSUj) gives rise to the same J q dq divergence 
as the curvature term for each F^""^ . Let us evaluate its coefficient. Going 
from one side of the surface S to the other changes the dihedral angle 6 into 
271 — 9. The sum of the wedge contributions to 5p from the two neighbouring 
domains is thus 



127r2 



ds 



(vr -e){Ti- 59) {9 - vr) {59 - 97i) 
9 ^ 271-9 



1 f.A^~(^) 



2 



— ds 



6n J 9(271-9) ' 
(80) 

and it is a positive number as soon as 6' 7^ tt. Hence Freg is dominated by 
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a positive, divergent term proportional to (jSUI) . and the field exerts in this 
idealized model an infinite constraint on the considered foil, which tends to 
flatten its dihedron. 

The same divergence occurs for an open conducting foil. Its edge is 
equivalent to a dihedral angle 9 = 27r, the contribution of which to p^^^ 
is (3/87r) J ds. The divergence of F implies, for instance, that the zero-point 
energy of the field produces a strong attractive force which tends to join the 
two halves of a thin conducting foil that is cut along a line. 

We also find a divergence if S involves several adjacent dihedra, for in- 
stance, if it is made of three half-planes joined along their common edge as 
in a honeycomb. In this case, for 9 = 2tt/3, we find a contribution to 6p 
equal to (— 7/247r) f ds which is negative. Such a configuration would thus 
be particularly stable if the experiment could be realized. 

The fact that one cannot define a finite Casimir energy for thin conducting 
foils with creases is related to the singularities of electromagnetic fields near 
sharp edges. 

Nevertheless the above treatment can easily be adapted to Casimir forces 
between perfectly conducting rigid bodies, even if they have sharp angles. 
Consider, for instance, two bulky wedges. The free energy of each of them 
cannot be renormalized, because the kernel K{a, j3) is singular as l/\a — /3p 
for two neighbouring points a and (3 located on the different sides of the edge; 
the resulting contribution TiK^ to (j78|) . in particular, is seen to generate a 
divergence in (|76|) . If, however, according to the remark at the end of section 
15.51 we focus on the interaction between the two indeformable wedges, letting 
aside their own energies, such divergent contributions are irrelevant. The 
interaction free energy of the two wedges is thus obtained by keeping, in the 
expansion of (j78j) in powers of K, only those terms which involve scatterings 
on both wedges (with an even number of factors K). Their contribution 
to (fTBj) is expected to be finite in spite of the short-distance singularity of 
K near the edges of the wedges, and the Casimir force between them is 
therefore well defined. For example, for two wedges with the same dihedral 
angle 29 facing each other perpendicularly at a distance L, the two-scattering 
approximation expressed below by eq.(|92|) provides at T = an attractive 
Casimir interaction energy equal to —he ig^9 / 'in'^L. 

We now review some applications of the general expression (|7B|) for the 
Casimir free energy. More details can be found in |18j . 
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6 Applications 



6.1 Parallel plates 

For two parallel plates with area A, lying at a distance L from each other, 
the function \E'(?/) associated with the three regions separated by these plates 
is 

ny) = ^Hi-^~'''^), (81) 

wherefrom we get the elementary Casimir effect at T = 0. The low temper- 
ature expansion of the free energy (f76|l . 

agrees with the result (jTFjl of the direct calculation, and shows that the 
Casimir attraction —OF/dL increases with the temperature. 

The comparison of (jZOI), which includes the factor ln(l — e~^'^/'^), with 
(irnjl . (jHH), exhibits a duality between high and low temperatures, as antici- 
pated in Eqs dni, ^ and (jni. 



f(T)-f(0)^(H^)^[f(^)-F(0) 



At high temperature, we have 



F{T) = - ^l^f^^ + OiT^-^^"^^'^) , 



(83) 



(84) 



yielding again an attraction due to radiation pressure. 

The entropy —dF/dT rises at low temperatures as 3^T^(^(3)/27r/),^c^, 
independently of the distance between the plates, and tends to a finite limit 
ACii) I^T^L"^ at high temperature. 



6.2 Low temperatures 

Owing to the more and more rapid oscillations of g{y) when T ^ 0, we can 
evaluate the second term of (f7H|l for low temperatures by expanding \E'(?/) 
around ?/ = 0. This yields 

^^^^ ' = 3^^^'^°^ - 135^^'"^°) + • ^^^^ 
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This behaviour is related to the topology of the boundaries S, since ac- 
cording to (f7^ \l/'(0) vanishes for a singly connected surface, and equals nA 
for a multiply connected surface. Accordingly, the low-temperature entropy 
arising from (jH3j) behaves for n = as (like the entropy of the black body), 
but is large as — 27r n AT/3hc for n ^ 0. This negative sign looks paradoxical. 
It is related to the fact that for torus-like topologies, permanent supercurrents 
can generate static magnetic fields. The occurence of a number 2n of such 
modes with q = 0, which do not contribute to the Casimir effect, entails a 
depletion in the distribution p^'"\q) for g 7^ 0. In fact, we see from (f7T|). 
(f75j) that 5p{q) for S — > 00 tends to — 2nA/7r as g — 0; it is this negative sign 
which is reflected in that of the Casimir entropy. However the total entropy 
of a quantum system must be positive. In the present case this results from 
the positivity of the total density of eigenmodes p^'"^ = p^'^^ + 5p of the 
field, ensured by the fact that Vq^/ir'^ is infinite in the large S limit 

considered here. 

The dimensionless parameter of the expansion ()85|) is IT /he where I is 
the typical size of the system S. The lowest order contributions should 
become experimentally accessible since this parameter is 0.5 for T = 300K 
and £ = 3/im. 

The different behaviour of ()82|) and ()85p arises from the fact that ()85p 
holds only for a finite system, whereas for an infinite system like a pair of 
parallel plates ^'"{y) diverges when ?/ — as shown by (jHTj) . 

6.3 High temperatures 

At high temperature the second term of ()76p is dominated by the first saw- 
tooth of giy). The corresponding calculation yields 

F = -CT\n{T/hcQ) + 0{T-^) , (86) 

c 5^ //"(F-sk)-"' f'^' 

lnQ = -- dy\ny^>'{y) . (88) 

Here the dimensionless parameter hc/TR of the expansion ()86|) is governed 
by a length R of the order associated with the short-range behaviour 
of the kernel K. This characteristic length is therefore a typical curvature 
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radius R of S. The high-temperature hmit might become experimentally 
relevant for crippled foils with small R. 

The dominant term of ()86|) is formally the same as the free energy of 
a number C of classical harmonic oscillators with average frequency cQ/27r. 
Contrary to what happens for the black-body which requires a quantum 
treatment at any temperature, the Casimir contribution that we have calcu- 
lated, which describes the change in the free energy of photons brought in 
by the boundaries S, takes a classical form for T S> he/ R. This is possible 
here because the modes that contribute to the Casimir effect have bounded 
frequencies, whereas the modes with hv ^ T crucially contribute to the 
black-body radiation. 

The internal energy U ~ CT arising from (|86p expresses the classical 
equipartition, and as usual in the classical limit the entropy Cln(eT/^cQ) 
depends on Planck's constant in its additive constant. The number C, pos- 
itive or negative, is interpreted as the average number of modes with finite 
frequency added by the introduction of the boundaries S. Here again the 
topology of S enters the expression (jHTj) of C through the genus n and the 
integer J c?^a/47ri?ii?2- For parallel plates, we have C = but ClnQ = 
C(3)/87rL2. 

The high-temperature Casimir constraints on the conductors S describe 
the effects of radiation pressure. To dominant order they behave as TlnT 
and are obtained by studying how C varies when the conductors are displaced 
or deformed. They tend to let F decrease, thus to let C increase. The only 
non-topologic part of C, 3/ d'^a/{32nR'^), does not depend on the relative 
position of the conductors. Hence, to the dominant order in TlnT, there 
are no forces between different conductors induced by the field. Moreover, 
since C is dimensionless and scale-invariant, there are no there forces tending 
to dilate or contract hollow conducting shells. However C increases with 
the average curvature of S, so that the Casimir effect tends at high 
temperature to let conducting foils undulate. This tendency is limited by the 
next terms of the expansion (jHBj) . and the curvature |-R|~^ tends to rise up 
to values of order T/hc. 

The next order contributions to the constraints, of order T, arise from 
the contribution CT In Q to F. Since is proportional to the size of S, 
these Casimir forces of order T would tend to contract S* if C is negative, to 
expand it if C is positive. 
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6.4 The wrinkling effect 

The existence of constraints that tend to wrinkle conducting surfaces at high 
temperature is confirmed by the study of small deformations of a thin foil. 
The evaluation of (fTBj) for a weakly deformed conducting plane S is achieved 
by means of a two-dimensional Fourier analysis. It is found that the con- 
straints created by the field tend to create ripples with wavelengths larger 
than 2.9hc/T, and to restore flatness for smaller wavelengths. In particular, 
the Casimir effect proper, at zero temperature, tends to suppress the curva- 
tures. Thus a conducting plane foil is stable atT = 0, but unstable atT^O 
under small deformations. 

This phenomena is confirmed by the study of the distribution in space of 
the free energy of the electromagnetic field. The density of Casimir free energy 
/(r) is obtained in the same way as the total free energy (fTBj) . except for the 
integration over the point (jUH|) . In fact, taking into acocunt the spectral 
representation ()61|) of the magnetic Green function F and the similar one for 
the electric function S, the contribution of each pole k = oi ^^"\k) is 
then weighted by i[B^(r) + E^(r)]. (The coefficients /ig ^ and eo of are 
recovered when the physical fields entering ()45|) are expressed in terms of the 
real and normalized functions entering T and H.) The density of free energy 
/(r) is thus found as a series describing a wave starting from r, scattering 
an even number of times on S and returning to r. The two-scattering term 
is sufficient to provide the free energy density at a point r located near tS*, at 
a distance d from it much shorter than the local curvature radii. We find at 
low temperature 

tic T^Ci^) 



and at high temperature 

T 



fir) 



IQnRd'^ 



\n{2dT/hc) + C-^ 



(90) 



where C is Euler's constant. 

This density is not bounded when d ^ 0, and its divergence is non- 
integrable. Hence, even after renormalization by subtraction of the free en- 
ergy of the vacuum without boundaries, the free energy F^""^ associated with 
each region is divergent. The total free energy F = ^^^^ finite because 
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on the two sides of S the average curvatures are opposite at each point, 
so that the divergences from (|89|) or from (|9U|) cancel each other. 

The sign in ()89|) shows that the presence of a perfectly conduting foil 
produces the transfer of an infinite amount of zero-point energy from the 
concave to the convex side, whereas, according to (jHOJ, the energy of real 
photons is transferred from the convex to the concave side. These opposite 
signs are consistent with the stability or unstability of a plane foil against 
deformations, depending on the temperature. 

6.5 Other examples 

6.5.1 Van der Waals and Casimir- Polder forces. 

The Casimir forces between conductors lying far apart can be evaluated by 
means of the free energy (fTBj) of the vacuum separating them. We noted 
that it may also be attributed to the random currents that circulate on their 
surface and produce the field. Such forces have thus the same nature as van 
der Waals forces of mutual induction, except for the retarded character of 
the interaction. We find from (fTBj) that two conductors at a large distance L 
apart attract each other as at zero temperature (Casimir- Polder forces 
jH]), as T/L'^ at high temperature (van der Waals forces). Torques are also 
found for anisotropic bodies. The same results hold for a small conducting 
body facing a mirror, which is attracted by it |31j. 

6.5.2 Derjaguin approximation. 

For two neighbouring conductors, the multiple scattering expansion of (f75|l 
can be used to justify the Derjaguin approximation seen above. Consider 
for instance the Casimir force between a plane and a sphere with radius R, 
and denote by L their shortest distance (Fig. [7j). The trace in (f75|l can be 
replaced by an integration over a point a of the plane: 

^{y) = -\ld^ay^[trHl-K^)l^ , (91) 

where the trace tr is meant on the tensor index only. Suppose L <^ R. 
Owing to the exponential decrease of K, the integral is dominated by 
the contributions such that a lies at a distance x from the sphere of order L, 
and such that all successive scattering points also lie at a distance of order 
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Figure 7: Multiple scattering contributing to the characteristic function \& 
(ED). 

L from a. Thus, for each a, the integrand of ()91|) is approximately the same 
as for two parallel plates lying at a distance x apart, which according to 
fl81|) is given by 7r~^?/2ln(l — e"^^^). This is just Derjaguin approximation. 
Corrections can be obtained from (|9ip. 

6.5.3 Two-scattering approximation. 

Another useful approximation is the two-scattering approximation^ for which 
we retain for '^{y) only the lowest order term ^TrKydK/dy of (fTSj) . More 
explicitly, this yields (Fig. ^ 

¥^\y) = 2y^ [ d2^d2^dGo{\a-P\)dGo{\a-f3\) 
dy J dua dup 

= -— / d'^ad'^l3{n^.p){np.p)-- — 

Stt^ J p dp 

where p is the vector a — (3. Numerical tests show that this approximation 
should be fairly good; for example it yields for the Casimir force at T = 
between two parallel plates the correct result times QO/vr^, an error of 8%. 

Using this approximation, we can evaluate the free energy of interaction 
between a plane and a sphere. The integrals in (jU^ can be completely worked 
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Figure 8: Geometry of the two-scattering approximation ()92j] . 



out and yield 



= + 1) e-'y^ - (^-Ry + e-'y^-'y^ . (93) 

Hence we find at T = 0, for L -C i?, 

E^^' ^ H , (94) 

which exhibits a correction in L/i? to the two-scattering Derjaguin contribu- 
tion. 

The study of a spherical shell S with radius R shows that the Casimir 
energy at T = behaves as HE 

E = 0.046 hc/R . (95) 

At high temperature, we find in agreement with ()86|) 



F=--|ln(™/M+ 0.769] -^-j _ + O j . (96) 

In the previous examples we had found attractive Casimir forces. In contrast, 
these forces tend here at any temperature to expand the sphere. They increase 
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with T. The radiation pressure exerted from inside thus exceeds that exerted 
from outside, contrary to what happens for parallel plates. 

An intermediate geometry between a sphere and parallel plates is that of 
a cylinder. We evaluate the Casimir energy of a hollow cylinder S by means 
of the two-scattering approximation (jU^ . For a cylinder with radius R and 
length C ^ R, we find at low temperature 

F^'^ ~ -TfS^^ , (97) 

45 h^c^ ^ ' 

and at high temperature 

in agreement with (jHE))- The cylinder tends to shrink at high temperature. At 
zero temperature, the zero-point energy vanishes in the considered approxi- 
mation and its exact value is very small, an intermediate situation between 
the parallel plates (|H^ and the sphere (|^. 

Unfortunately all these Casimir constraints on thin conducting foils, in- 
cluding the wrinkling effect, are presently difficult to detect experimentally, 
because of their weakness compared to the binding forces that ensure the 
cohesion of the metallic sheets. 
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